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Abstract 

We review the phenomenology of electro-weak bosons produced in hadron-hadron 
collisions. In particular, we discuss the transverse momentum distribution of lepton 
pairs with invariant mass close to the Z peak. We describe the theoretical calculation 
for the magnitude of the transverse momentum Qt and its comparison to Tevatron 
and LHC data. We also discuss the related variable (f>* , describing its experimental 
advantages as well as its relation to the standard Qt variable. Finally, we compare 
resummed predictions for (jf to data. 
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1 Introduction 



The production of a lepton pair l + l~ in hadron-hadron collisions is one of the most extensively 
studied processes in particle phenomenology, with the original paper appearing more than forty 
years ago pp. Since then a huge theoretical effort has gone into improving the accuracy of the 
predictions. For instance, QCD corrections are known to next-to-next-to leading order (NNLO), 
not only for the inclusive rate [2], but also for rapidity distributions [3] and more recently for 
fully exclusive cross-sections [1HS]- 

Of particular interest is the transverse momentum distribution of the lepton pair, or equiva- 
lently of the gauge boson decaying into it. At Born level, the production process is described by 
quark-antiquark annihilation diagrams. In the framework of collinear factorisation, the initial- 
state (anti) quarks are on-shell and they move in same direction as the incoming hadrons. Thus, 
at this perturbative order the vector boson (Z/7*), and hence the lepton pair, is produced with no 
transverse momentum with respect to the beam. The vector boson acquires a non-zero transverse 
momentum when we consider QCD radiative corrections. The Z/'y* transverse momentum distri- 
bution is therefore sensitive to multi-gluon emission from the initial-state partons, while having 
a particularly simple final state, and hence provides a powerful tool to test our understanding of 
QCD dynamics. 

The observable we are consideing is a typical example of a multi-scale problem and the correct 
treatment of multi-gluon emissions goes beyond fixed-order perturbation theory. Let us introduce 
the invariant mass of the leptons M , that will be chosen around the Z mass, and the let us call 
Qt the magnitude of the vector boson transverse momentum. We are going to consider three 
different regimes. When Qt ~ M we expect fixed order perturbation theory to work. In the 
region Aqcd "C Qt <S M, we can still rely on perturbation theory but large logarithms of the 
ratio Qt/M may spoil the convergence of the perturbative expansion and must be resummed to 
all orders [9H19], Finally, in the region Qt ~ Aqcd we expect non-perturbative effects to play 
a significant role. These contributions are associated with the intrinsic transverse momentum of 
the initial state partons, which is neglected in the usual perturbative treatment based on collinear 
factorisation. Therefore, it is important to compute a solid perturbative prediction, so that we 
can compare it to precise data coming from the experiments and be able to eventually pin down 
non-perturbative contributions. 

2 The challenge of theoretical predictions at small Qt 

The leading order (LO) transverse momentum distribution can be obtained by integrating the 
0{a s ) corrections to the Born process [8]. In the small-Qr limit we obtain the following behaviour 



The distribution above exhibits a double logarithmic divergence when Qt — > 0, which corresponds 
to the emission of a soft and collinear gluon. The coefficient of the double logarithm (DL) is very 
simple: A = 2C^Xo, where So is the hadronic Born cross-section. The single-logarithmic term 
instead contains a contribution which originates from the emission of a hard collinear parton. For 
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this reason the coefficient B does not trivially factorise and it contains a convolution term which 
depends on the parton distribution functions (PDFs). 

The behaviour described at 0(a s ) by Eq. (pQ) is actually present to all orders in perturbation 
theory. At the DL level we have 



da 



dQ T Q2 



(2) 

rp / ^T **£*H ' ^T 



Not only is any truncation of the above series divergent in the formal limit Qt —> but it 
also loses predictive power if a s In jL- = 0(1), because all higher-orders are equally important. 
These contributions need to be resummed to all-orders in order to obtain a reliable answer. The 
calculation in the soft and collinear limit is not difficult and we find that the coefficients Ai are 
proportional to powers of the Casimir Cp. The DL series can be then summed to: 

da d -if^in 2 ^ 

The resummed expression exhibits a Sudakov form factor that suppresses the emission of soft 
gluons. As a consequence, the distribution is now well-behaved and, actually, vanishes in the 
Qt — > limit. However, we may question the usefulness of Eq. Q for phenomenology. In 
fact, the hypotheses used to derive the DL results are too crude to capture all the relevant 
effects, though they are formally sub-leading. For instance, using momentum conservation in the 
transverse plane, we have Q T = — X^Ati) where are the transverse momenta of the emitted 

partons. It is then clear that small values of Qt = \J\Qt\ 2 can ^ e acn i eve ol via two distinct 
mechanisms: emission of gluons with small transverse momentum or kinematical cancellation 
between, not necessarily soft, k_ T i- 

We would then like to extend the resummation formalism beyond the DL level. The derivation 
we sketch in the following is valid at next-to- leading logarithmic accuracy (NLL), although the 
state of the art is actually one logarithmic order higher (NNLL). Firstly, we use the well-known 
property that convolutions of matrix-elements with PDFs can be written as ordinary products of 
their Mellin moments. Therefore, we are able to arrive at the distribution for the emission of an 
arbitrary number of collinear and optionally soft gluons in Mellin space: 

1 n=0 i=l ' \i=l / 

where P qq is the Altarelli-Parisi splitting function and £o(A0 is the Mellin transform of the 
hadronic Born cross-section 

Eo(JV) = d ^/ gM (iV,/i)/ g - /B (iV,/x). (5) 

We note that the 5 function that enforces momentum conservation in Eq. spoils the 
factorisation properties of the result. This issue can be solved by introducing a two-dimensional 



2 



impact-parameter representation of the 5 function 



« 2 > (i>,+2r)=i/ d 'v^n 

\i=l / J 8=1 



(6) 



which has the desired factorised form. We can now perform the sum in Eq. obtaining an 
exponential: 



dd 



dQ 



^oW^ / d V*P 



(7) 



where the (—1) term has been included to account for virtual corrections. If we are interested in 
the Qt spectrum, we integrate over the angle between b and Q T , obtaining a Bessel function: 



goCjV) 
4vr 2 



dbb 



2k 



dip e 



ibQ T cos ip -R N (b,M,N) 



£o(AQ 
2tt 



dbb,J (bQ T )e 



-R N (b,M,N) 



(8) 



We can now have a closer look at the resummed exponent Rn; at NLL accuracy we find the 
following structure: 



R N (b,M,N) 



M 2 



dkj. a s (kT) 



-a 7T 



7T J Kj, 



(9) 



where 6 = be YE /2 and we have introduced the coefficients: 



A( 1 ) = C F , A( 2 ) = ^ 



HZ _ E 

18 6 
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71/ 



(10) 



and 7 gg is the Mellin transform of the LO splitting function. The N dependent part of the function 
Rn can be absorbed into the PDFs using DGLAP evolution. The PDFs, which in Eq. (j3J) are 
contained in the function So, are then evolved down to the scale fi = 1/b. In the following the 
part of Rjy which does not depend on N will be denoted by R. 

At NLL accuracy the integration over fey in Eq. ([9]) must be performed with the two-loop 
expression for the running coupling. We can then express the resummed exponent as the sum of 
two contributions: 

R(b) = LgV(a s L) + gW(a s L) + ... (11) 

with 



5 (1) (A) 



AM 



7T/3 A 

5 (2) (A) =-^ln(l-A) + 
irp 



A-ln(l-A)] , 

A( 2 )[A + (1- A)ln(l-A)] A^p l 



(12) 



A) 



7T/3 3 



A±ML^) + Ii n2( i_ A) 



(13) 
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and A = a s (M)(3oL, L = In (6M) 2 . The function g^ resums LL terms, while NLL ones. 

As mentioned before, the state of the art in Qt resummation is actually NNLL. In order to 
achieve that accuracy the coefficients and [20] have been calculated. Moreover, the 

integrals in Eq. @ must be performed with the three-loop expression for the running coupling. 
This results in a new contribution, a s g^ , to Eq. (|lip . 

So far we have dealt with all the logarithmic terms in the function R. However, there is 
still one contribution we must consider in order to achieve (N)NLL accuracy: non-logarithmic 
corrections to the Born cross-section Eq in Eq. ([5]): 



1 a,P 



1=1 V / 



+ £ 

i=i 



7T 



f a/A (N,l/b)fp /B (N,l/b), 



(14) 



in particular, the i = 1, 2 coefficients are needed pT9|[2Tj . 

Lastly, in order to obtain a resummed expression for the Qt distribution, the integral over 
the impact parameter b in Eq. (J7J) must be performed. A closer look at Eq. (|lip reveals that this 
integration is problematic both at small b and large b, because the function R diverges in both 
limits. The small-6 region corresponds to large values of Qt, which are beyond the jurisdiction of 
the resummation formula. Different prescriptions that preserve the logarithmic accuracy can be 
found in the literature. For instance, in Refs [18^119] the argument of the logarithms is modified: 
In (b 2 M 2 ) ->■ In (l + b 2 M 2 ), so that R(b) vanishes as b ->• 0. Alternatively [22H24] we can freeze 
the radiator below some value 6 m j n . 

The large-6 divergence, a s /3oL = 1 is due to the Landau pole in the running coupling and is 
associated to non-perturbative (NP) behaviour. In this case as well different prescriptions to deal 
with the singularity can be found in the literature. In Ref. [12] a smooth transition between a 
perturbative and NP region was derived in terms of functions that must be fitted from the data. 
Improvements of this NP form factor were suggested in Refs. [25|. l26j . Other approaches consist 
in introducing an upper limit to the b integration [22H24] or deforming the contour integration 
off the real axis, so to avoid the Landau singularity, as, for instance, in pT8|[T9] . 

Different, and equally successful, approaches to Qt resummation are the ones based on soft- 
collinear effective theory (SCET). Calculations have been performed to NNLL accuracy by two 
groups in [20l[29] and in [30H32] and successfully compared to Tevatron and LHC data. 



3 Comparison between theory and experiment 

The transverse momentum spectrum of vector bosons produced via the Drell-Yan mechanism has 
been extensively studied by the Tevatron experiments [33H35] and, more recently, by the LHC 
collaborations as well |86,37J. State-of-the-art-perturbative predictions have been compared to 
data in order to test the accuracy of resummed calculations and to asses the role of NP effects. 

The phenomenological tool mostly used by the experimental collaborations is Resbos [25|.I26| . 
This program uses the resummation formalism developed in Ref. |12] in order to produce NNLL 
resummed distributions. Moreover, it is matched to fixed-order calculations, in order to provide 
reliable QCD predictions for all values of Qt- NP effects are also modelled in this code essentially 
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Figure 1: Comparison of the theoretical prediction from Resbos for the Qt spectrum to the 
experimental data collected by the D0 collaboration [35J, for the inclusive sample (on the left) 
and for forward rapidities (on the right). This figure is taken from Ref |35j . 



with a Gaussian smearing (BLNY form factor [26]), with parameters fitted to the experimental 
data. Additionally, phenomenological studies in semi-inclusive deep-inelastic scattering processes 
suggest an effective dependence of the NP effects on Bjorken x. This can be then translated 
into a small- x contribution to the NP form factor in Drell-Yan processes [27], usually referred 
to as small- x broadening. In Fig. [TJ we report the comparison between the Resbos theoretical 
predictions and the data collected by the D0 collaboration at Tevatron Run II |35j. The red 
line is the theory prediction obtained with the BLNY form factor in Resbos, for the inclusive 
sample, on the left, and for forward Z boson rapidities, on the right. In the latter case, also the 
prediction with small- x broadening is plotted. The agreement between theory and experiment 
is remarkable. Unfortunately, even at large rapidities, where small- x effects are expected to be 
more pronounced, we cannot discriminate between the two different NP models because of the 
experimental uncertainty. 

A complementary analysis of the D0 data was carried out in [19] and the results are shown 
in Fig. [2j The theory prediction is a pure perturbative calculation which resums logarithms at 
NNLL accuracy and it is matched to a NLO calculation. The band represents an estimate of the 
perturbative uncertainty of the theoretical calculation and, within this error, a good description 
of the data is found, without any NP correction. The theoretical uncertainty is computed by 
varying the arbitrary scales present in the calculation. As usual, factorisation and resummation 
scale variations are considered; however, it turns out that the dominant source of uncertainty 
at small-Qr comes from missing higher logarithmic terms, which are estimated by varying the 
argument of the logarithms we are resumming. A pure NNLL+NLO, without any explicit NP 
effect also gives a good description of LHC data [24], as shown in Fig. [3] 
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Figure 2: Comparison of the theoretical NNLL+NLO prediction of Ref. |19j for the Qt spectrum 
to the experimental data collected by the D0 collaboration. This figure is taken from Ref. |19j . 




Figure 3: Comparison of the theoretical NNLL+NLO prediction [23] for the Qt spectrum to 
the experimental data collected by the ATLAS collaboration [36J (on the left) and by the CMS 
collaboration [37] (on the right). This figure is taken from Ref. |24| . 
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4 Novel variables 



In the previous section we have seen that, despite very accurate calculations, no clear conclusions 
can be drawn by the comparison of the theory to the Qt data from the Tevatron experiments. On 
the theoretical side, we have different ways of estimating the uncertainty as well as different NP 
models. On the experimental side, we would like to reduce the uncertainty associated with the 
measurement. To this purpose new variables, labelled ot and (jf , were introduced in Refs [381139]; 
the <f>* distribution was measured by the DO collaboration [40] . 



4.1 Experimental motivations 

The variable ay is component of Qt orthogonal to the dilepton thrust axis [38j : 

a T = -=> p— , (15) 

\lTX — lT2\ 

where all the above three-vectors are defined by v = (v, 0) and l T i 2 are the leptons' transverse 
momenta. One of the main experimental uncertainties comes the resolution in the measurement 
of these momenta. We note that at low Qt we have: 

a T ^ -. —r— sin Acf>, (16) 

Iti + iT2 

where A<f> is the dilepton azimuthal separation, which is close to ir in the limit we are considering. 
Thus, at \ow-Qt, the uncertainty for ay is reduced with respect to the one for It% (and hence Qt) 
by the presence of the small factor sinA^ [38] . Moreover, it was noticed that this uncertainty 
can be further reduced if we rescale ot by the dilepton invariant mass M , essentially because 
resolution effects partly cancel in the ratio Ideally, we would like to define a variable that 
depends only on very well-measured angles and, at the same time, maintains a close relation to 
Qt- The observable (ft* satisfies these properties 



tan I IL_^f)sin0*, (17) 



where 8* is the scattering angle of the dileptons with respect to the beam, in the boosted frame 
where the leptons are aligned. This angle can be related to the pseudo-rapidities of the leptons: 
sm8* = 1/ cosh (7/' 1 ) — ?/' 2 )), hence (jf is fully determined by angular measurements. Moreover, 
in the Iow-Qt limit (jf reduces precisely to clt/M. 



4.2 Phenomenology of (jf 

We would like to derive a resummed expression for the (jf distribution and compare it to the 
DO data |40j . We have already seen that, in the relevant Qt — > limit, this variable reduces to 
cit/M, i.e. one component of the transverse momentum vector. The resummation for ar was 
discussed in [28] and it can be related to the traditional Qt resummation described in section [2j 
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However, in the present case, we are only interested in one component of Q T , rather than its 
magnitude. This difference results into a cosine rather than the Bessel function Jo in Eq. ([8]) 

J -i poo 

i=- = -/ dbM cos(bM<j)*)e- R ^ M ^(x u X2,b,M) . (18) 
d^* vr J 

We first note that the function R which resums the large logarithms in &-space is the same as 
the one obtained for Qt resummation Eq. (fTT]) . Secondly, the presence of the cosine function has 
important phenomenological consequences. In fact, the <j)* distribution does not show the typical 
Sudakov behaviour as Qt, but rather tends to a constant plateau as <p* — > 0. We have already 
discussed that low values of Qt (and hence also (jf) can be obtained via Sudakov suppression or 
kinematical cancellation. In the present case of (jf the kinematical cancellation starts to dominate 
prior to the formation of the Sudakov peak, in contrast to the Qt case. 

The D0 collaboration compared their data to theoretical predictions obtained from the pro- 
gram Resbos. Thanks to the smaller experimental uncertainty the collaboration was able to 
discriminate between two different NP models, showing in particular that the data at forward 
rapidities disfavoured small- x broadening, which had not been previously possible due to errors 
on the Qt spectrum, even with Tevatron Run-II data |40J . An independent calculation for the 
4>* distribution was performed in [23J. NNLL+NLO perturbative predictions were compared to 
D0 data with a good agreement in all rapidity regions, once the theoretical uncertainties were 
faithfully estimated. As an example, we report the comparison between data and theory in Fig. H] 
in the case of electrons, for central (left) and forward (right) rapidities. 

Predictions for the (jf distribution in proton-proton collision at 7 TeV were first computed 
in |24| . where measurements by the LHC collaborations were encouraged. 

5 Conclusions 

We have reviewed recent phenomenological studies of the transverse momentum distributions of 
vector bosons produced in hadron-hadron collisions. We have started by discussing the theoretical 
issues we need to face in order to produce a reliable prediction. The Qt spectrum can be computed 
in perturbation theory, but it is affected by large logarithms in the small-Qr region. These terms 
can be resummed to all-orders and when the result is matched to fixed-order calculation, we can 
obtain reliable predictions for all values of Qt- 

We have reported on the comparison between state-of-the-art theoretical calculations for 
the Qt spectrum and data coming from the Tevatron and the LHC experiments. In order 
to reduce the experimental uncertainty, the variable <fi* was introduced and measured by the 
D0 collaboration. We have discussed its relation with Qt and, consequently, its resummation. 
The very low experimental systematics makes the variable (ff invaluable for precise studies of 
electro-weak bosons at LHC, with the possibility of probing many subtle QCD effects, such as 
multi-gluon dynamics and its description within, or beyond, perturbation theory. 
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Figure 4: Comparison of the theoretical NNLL+NLO prediction [23] for the (f>* distribution to the 
experimental data collected by the D0 collaboration |40j . in different rapidity bins. This figure 
is taken from Ref. |23j . 
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